Rotating Morse wave packet dynamics of diatomic molecule 
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We investigate the dynamics of a rotating Morse wave packet, appropriate for a ro-vibrating 
diatomic molecule. The coupling between vibrational and rotational degrees of freedom is explicated 
in real position space as well as in phase space Wigner distribution of a SU(2) coherent state at 
various dynamically evolved times. We choose the well studied I2 molecule with the parameter 
values in good agrement with experiments. A quantitative measure of the angles of rotations for 
different angular momenta is also given. 
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PACS numbers: 03.65.Ge, 42.50.Md, 34.20.Cf 



I. INTRODUCTION 



The Morse potential [l[ is extensively studied in the 
literature and has been a subject of considerable atten- 
tion, ft involves the vibrational degrees of freedom of 
a diatomic molecule. On the other hand, rigid rotor 
model is described to take into account only the rota- 
tional degrees of freedom, where on-axis rotation is usu- 
ally neglected and end-over-end rotations are identical for 
a homo-nuclear molecule. The rotational and vibrational 
energy scales are far apart and it has been a customary 
to treat them independently. However, in reality, this 
separation is inadequate and the internuclear distance is 
affected by the centrifugal force resulting from rotation. 
Thus, one should consider both degrees of freedom of a 
diatomic molecule. A combined form of the Morse po- 
tential and the rotational centrifugal term is known as 
the rotating Morse potential. It leads a more complex 
energy spectrum reflecting the coupling among the two 
intrinsic properties of nuclear molecule. It is known that 
the Schrodinger equation for this potential can be ex- 
actly solved only when the orbital angular momentum 
quantum number j is equal to zero. Since, there is no 
exact analytical solution for the rotating Morse potential 
in general, some semiclassical and/or numerical solutions 
have been obtained by using various approximate meth- 
ods. The original approach developed by Pekeris Q is 
based on the expansion of the centrifugal barrier in a 
series of exponentials depending on the internuclear dis- 
tance up to second order. A more appropriate expansion 
has been proposed by Duff and Rabitz [U and Elsum 
and Gordon . 4| . These semi-analytical models are based 
on the idea of building up Morse potentials with parame- 
ters dependent on the rotational quantum number, which 
represent adequately the effective potential curves of the 
rotating Morse oscillator. There are various other numer- 
ical and proposed methods to study the rotating Morse 
potential @, i, B 1, [1 EI El El El ■ These studies are 
very useful to obtain more accurate energy eigen values. 
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Recent development of ultrashort pulses opened a new 
research field which allows to study the wave packet dy- 
namics of diatomic molecule in real time scales. Pulses of 
a few tens of femtosecond allow a coherent broadband ex- 
citation with preparation and detection of rovibrational 
dynamics. In all previous works, the wave packet dynam- 
ics have been studied in diatomic molecule either consid- 
ering the vibration 0, EE Ell and rotation [H El, El 
separately. There are few works, which considers the 
both in two-time scale configuration [H, [U [11] . 

In this work, we study the wave packet dynamics, con- 
sidering the coupling among the vibrational and rota- 
tional degrees of freedom of nuclear motion. Our an- 
alytical approach for the study of the rotating Morse 
potential is inspired by Gordon's method Q, where the 
effective potential becomes minimum around a new equi- 
librium internuclear distance rj, which is a function of 
angular momentum quantum number j. This model, 
so far, most appropriately manifests the coupling be- 
tween two degrees of freedom, vibration and rotation. 
To the authors' knowledge, there is no study about the 
wave packet dynamics of this appropriate system. We 
show that this ID model is accurate upto a reasonably 
high angular momentum (j — 200), accessible to exper- 
iments. We study rovibrational interplay through a co- 
herent state (CS) dynamics, not only in coordinate space, 
but also in phase space Wigner distribution. It is shown 
that rovibrational interplay is better explicated in phase 
space rotation of the CS for different angular momenta. 
When rovibrational dynamics is seen in 3D configuration 
space, the corresponding phase space is six dimensional, 
which is not possible to visualize. Avoiding this 3D intri- 
cacy, this ID rotating Morse potential can well describe 
the coupling between the two degrees of freedom, hence 
can nicely manifest the rovibrational interplay for this 
realistic system in phase space. It also shows the ef- 
fects of rotation-vibration in the quantum interference 
regime. We have chosen the iodine molecule (I2) which 
is an unique suited seed molecule for laser induced flu- 
orescence. Additionally, we give a numerical method, 
which can estimate the angles of rotation, are in good 
agreement with the Wigner plots. 

The paper is organized as follows. In the following sec- 
tion, we present a brief overview and justification of the 
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FIG. 1: (a) The variation of rj with the rotational quantum number j. It implies that the numerical values, obtained by solving 
the transcendental Eq. [red filled circles in (a)] matches nicely with the approximate values from Eq. ((3| [black circles in 
(a)]. It starts to differ for higher values of j (> 200). The change of the dissociation energy of the effective potential with j is 
depicted in (b). 



ID rotating Morse system. We also evaluate the corre- 
sponding SU(2) CS wave packet. In Sec. Ill, we analyze 
the dynamics of the CS of a diatomic molecule and ex- 
plain the effect of rovibrational coupling in this. After 
showing the effect in the configuration space, we move 
to the phase space picture for a better representation, 
which clearly displays the rotational effect. A quantita- 
tive analysis of the angle of rotation is also given. This 
numerical analysis confirms the values of rotation angle, 
corresponding to various angular momentum quantum 
number. Finally, we end with some conclusions in Sec. 
IV. 



II. ROTATING MORSE WAVE PACKET 

The Schrodinger equation for the rotating Morse oscil- 
lator reads as follows, 



where A = ^y 1 ^ . Alternatively, rj can also be com- 
puted numerically by solving the transcendental equation 



"2^ + ^ (r) 



(1) 



where fi is the reduced mass of the diatomic molecule and 
n, j are the vibrational and rotational quantum numbers, 
respectively. The effective potential V e ff (r) has the form 

V eff (r) = D[e- 2 ^- r °) - 2e-^ r - r ^} + M + 1 } h \ ( 2 ) 

2nr 2 

The first part describes the well known Morse potential, 
where D is the dissociation energy of the Morse potential, 
r is the equilibrium internuclear separation and (3 is the 
range parameter. The second part of the above equation 
stands for the centrifugal contribution of rotation. More 
accurate description of the system can be obtained with 
a modified equilibrium internuclear distance Tj and a dis- 
sociation energy Dj 4] . Using a semi-analytical method 
12311 . one can find 
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These two sets of rj's are plotted in Fig. QJa), which 
shows very good agrement for j < 200. It shows that the 
rotational potency results into the increasing value of the 
equilibrium distance [Fig. [T{a)] and decreasing nature of 
dissociation energy [Fig. [ljb)]. Physically, in presence 
of the rotational centrifugal force, the two constituent 
atoms of a diatomic molecule tend to settle in a larger 
distance and they are more prone to dissociate, reducing 
the amount of energy required to make them indepen- 



dent. We define Aa 



and expand the centrifu- 



gal term of Eq.Q around r = Tj. Keeping terms upto 
second order, the effective potential becomes 
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where y^2X j e-^ r - r «^2\ j ue-'^ r ~ r ^ and u=e~^- ro \ 
The parameter, Xj — where the constants, cq, Ci, 

and C2 are dependent on quantum number j: 

c = 3Ajb] - SAjbj + Aj, 
ci = {ZAjb 2 - 2Ajbj + Aj + uD)/u, 
and c 2 = {3Ajtf - Ajbj + Aj + u 2 D)/u 2 , (6) 



with bj = (/3rj) _1 . Under this approximation, the 
Schrodinger equation |T]) is solved to yield the eigen func- 
tions of this rotating Morse system: 



i> n Ay) = N n ,je-y/ 2 y s L 2 n s (y), 
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where the variable y: 0<y<oo,nis the vibrational 
quantum number, L 2 f(y) stands for associate Laguerre 



polynomial and s — 



The number of 
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bound states in the system is [Xj — 1/2], where [p] de- 
notes the largest integer smaller than p. The parameter 



Xj is related as Xj = f*Xj. 



They satisfy the constraint condition s + 2n = 2Aj — 1. 



N n _j is the normalization constant: 



Nn,j = 



/3(2Xj -2n- l)r(n+ 1) 



1/2 



(8) 
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The rovibrational energy eigen values E n j turns out as 



E n>j = 2^(n + 1/2) - ^(n + 1/2) 2 + c - ^. (9) 
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r e = 5.03 a.u. to Tj = 5.0323 a.u. and Tj — 5.034 a.u. re- 
spectively) and dissociation energy diminishes. The rovi- 
bratational coupling changes the effective potential and 
hence energy levels are shifted upwards. The 10th en- 
ergy level, where the coherent state is excited, is shifted 
from E w = -0.04719 a.u. to E w = -0.046595 a.u. and 
Eio = —0.046108 a.u respectively. In our further study, 
in particular, we consider these effective potentials and 
the initial wave packet as a coherent state involving few 
lower levels of I2 molecule centered around n = 10. If we 
take only bound states of the potential, the dynamical 
symmetry group is SU(2). In this case, the displacement 
operator coherent state can be constructed and written 
as 



It is worth pointing out that in absence of rotation cq = 
0, C\ — c 2 = D, the system only describes a vibrating 
diatomic molecule. 
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III. DYNAMICS OF A ROTATING MORSE 
WAVE PACKET OF A DIATOMIC MOLECULE 

Ultrashort laser pulses can probe the molecular dy- 
namics in real time scale. Time-resolved measurements 
on the femtosecond scale have made it possible to know 
the evolution of a coherently prepared superposition of 
ro-vibrational states. Here, we consider a rovibrational 
wave packet of I2 molecule, which is a coherent state, de- 
pendent on particular rotational quantum number. Many 
theoretical and experimental investi gati ons have been 
carried out on this molecule [22], [24l. |25|] and has been 
a subject of many femtosecond studies, particularly in 
wave packet interferometer [2(1 The correspond- 

ing parameter values are (3 = 0.9849 a.u., reduced mass 
fi = 11.56 x 10 4 a.u., tq — 5.03 a.u., and D = 0.057 a.u.. 
In particular, Zewail and co-workers investigated rovi- 
brational wave packet dynamics in the well-characterized 
electronic B0+ state [24]|. Lohmuller et al. [22j dis- 
cussed about the pump-probe experiment of I2 in a room 
temperature and the detection of fractional revivals us- 
ing full-dimensional quantum wave packet. Following 
Lohmuller et at, we consider an initial rovibrational wave 
packet centered around 10th vibrational energy level and 
for fixed j = 60. This can be excited by using a 620 nm 
pump pulse, which is in good agreement with the exper- 
iments [3, HH and for detection purpose, probe pulses 
of 310 nm and 330 nm are used. These parameters allow 
the probe excitation mainly between the B0+ and /0+ 
states [29]. Under the laser polarizations magic angle 
conditions, it takes into account the vibrational as well 
as the rotational motions. 

The rotational effect on the effective potential of I2 
molecule is depicted in Fig. O Without rotation i.e., 
j = (solid line), it is nothing but the Morse poten- 
tial. Once rotation is considered, e.g., j = 60 (dashed 
line) or j — 81 (dotted line), the potential changes, i.e., 
the potential-minima slowly shifts towards right (from 
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FIG. 2: Effective potentials for j — (solid line), j = 60 
(dashed line) and j = 81 (dotted line) are depicted. The 
corresponding 10th energy levels are zoomed in the inset plot. 
The wave packets at t — 0.25 T rev are also shown in the inset 
as dark filled (j = 0), dashed {j = 60) and light filled (j = 81) 
plots, respectively. Changes of j value cause slight changes in 
the effective potential, which result in rotational effect on the 
wave packet. The two splitted wave packets in co-ordinate 
space come closer with increasing j and superpose each other 
for j — 81, producing interference ripples. The potential and 
the internuclear distance are in atomic units. 



Here, the group generators are j dependent and the 
weighting coefficient becomes 



dl = 



(n' — n)\ 



!r(2A j 



n!T(2A j - n') 



(11) 



It can be easily verified that for j = 0, the above coher- 
ent state (CS) is the same as reference [3(| ■ The SU(2) 
generators satisfy the commutator relations, [J+, J-] — 
2J , [Jq, J±] = ±J± and also satisfy 



[J+,J-]^n,j(.y) = 2joipn,j(y), 



(12) 
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FIG. 3: (Color online) Time evolution of Wigner function of the CS for j = at times (a) t = 0, (b) t = T rov /8, (c) t = T rcv /4, 
and (d) t = T rcv /2. The second row shows the corresponding plots for j — 60 [see (e)-(h)]. These structures (j — 60) reveal 
the different rovibrational coupling at different evolved times. Position and momentum variables, r and p, are in atomic units. 



at the level of wave function, where jo — n — Xj + 1/2. 
Jo is the projection operator of the angular momentum 
m: m = n — \j + 1/2. Although, the CS of a physi- 
cal system is never ideal, like that of harmonic oscillator, 
we find that our chosen CS is very well localized and re- 
mains coherent for a reasonable time duration. Presence 
of quadratic term in energy expression, leads some inter- 
esting phenomena, called fractional revivals, which occur 
at some specific instances between two full revivals (T rev ) 
[3ll [HJ • The short time evolution displays a classical pe- 
riodicity (T c i). The classical and revival time periods are 
respectively given by 

Tel = = 27rX \ and T mv = 2n\ 2 /c 2 . (13) 
2ci-c 2 Aj J 

For I2 molecule, these time scales are T c i = 0.156 pi- 
cosecond and T rev = 36.2 picosecond. At fractional re- 
vival times (rational fraction (r/s) of the revival time 
T mv ), the wave packet breaks into several number of sub- 
sidiary wave packets, where r and s are mutually prime 
integers. For even values of s the wave packet breaks 
into s/2 in number, otherwise breaks into s parts. In 
the course of time evolution, one obtains the Schrodinger 
cat like state at one fourth of the revival time. Four-way 
break up or the compass-like state emerges at one-eight 
of the revival time. Latter produces sub-Planck scale 
structures in phase space Wigner distribution. In the 
inset of Fig. [3J the 10th energy levels for different an- 
gular momenta are zoomed and the rotational coupling 
effect is shown at t = 0.25 T rev , when CS is splitted into 
two parts. Initially (j — 0), the two parts are situated 
at r — 4.7 a.u. and 5.58 a.u. (dark filled plot). After 
the introduction of rotational coupling (J — 60), they 
come close to each other, situated at r = 4.87 a.u. and 



5.48 a.u. respectively (dashed line). For a greater angu- 
lar momentum (J = 81), the position space probability 
structure looks completely different and shows oscillatory 
structure (light filled plot). The interpretation lies in the 
fact that both the splitted CSs oscillate inside the poten- 
tial well in back and forth motion. In the first quarter 
of the oscillation, they approach each other, while in the 
next quarter, they recede. At halfway of the oscillation, 
they reflect from the potential well with a phase change tt 
and become again counter-propagating. For the angular 
momentum j = 81, they overlap each other in the course 
of their oscillation and produce the oscillatory ripples, 
clearly visible in the inset of Fig. [2] and in Fig. 0Ja). 

Until now we have confined our study in position space 
explanation. This explanation does not reveal the full 
description of rotational coupling of the system. Hence, 
we enlarge our view towards the phase space description. 
Wigner function [33[ is very useful tool to visualize the 
quantum dynamics in phase space and is defined as 

1 r +oc 

W(r,p,t) - — / $*(r-r',*) 

x<5>(r + r' ) t)e- 2lpr '/ h dr' . (14) 

This representation can reveal interesting structures of 
the CS at different evolved times. As I2 molecules has 
large number of bound states (117), we choose the coher- 
ent state parameter a = 2.15, which involves few lower 
levels of the potential. In Fig. 02 first row shows the 
time evolution of the CS with j = 0. It represents the 
Morse vibrational dynamics at different fractional revival 
times. The initial CS is situated near the left poten- 
tial boundary and is well localized [Fig. [3] (a)]. Owing 
to the anharmonicity of the potential, the probability 
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FIG. 4: (Color online) (a) Probability density of rotating 
Morse wave packet for j = 81 at t = T rov /4. (b) The corre- 
sponding phase space Wigner distribution at the same time. 
Position and momentum variables, r and p, are in atomic 
units. 




FIG. 5: Variation of angle of rotation <f> of the wave packet 
with the angular momentum quantum number j. The dotted 
lines specify two particular rotations corresponding to two 
angular momenta. <f> is in the unit of radians. 



densities near the left boundary are stiller compare to 
those near right boundary. From the previous explana- 
tion, it is evident that the compass-like state appears at 
t = Trov/8 [Fig. [3] (b)], which shows the smallest inter- 
ference structures i.e., the well known sub-Planck scale 
structures [3(| [H, HH, HU . Origin of these smallest struc- 
tures is the superposition of the two diagonal interfer- 
ence ripples. The Schrodinger cat-state [37| appears at 
t = T'rcv/4 [Fig. Efc)], where the interference ripples ap- 
pear in the middle as alternate positive and negative val- 
ues, along the line joining the two component CSs. All 
these negative phase space interference structures are the 
strong signature of quantum behavior. At half revival 
time the CS revives near the right boundary of the po- 
tential well [Fig. [3] (d)]. Upon introducing the angular 
momentum contribution, we show the phase space struc- 
tures at the same times as above in Fig. [3Je)-(h). The 
initial CS is similar because it does not depend on rota- 
tion. One can observe the rotational effect on the other 
structures. For heavy molecules, the rotational effect is 
less for smaller values of j. We choose j = 60 as it is 



reasonable and also appropriate for experiment. Inter- 
estingly, the revival time T rev , does not change with j 
values, because the ratio A|/c2 remains same. Whereas, 
the classical period T c i alters due to the rovibrational 
coupling. The rotation is anti-clockwise, thereby making 
the angle of rotation positive in our convention. Fig.[3Jf) 
describes the rotational effect at the sub-Planck interfer- 
ence structures. It is established that these structures are 
very sensitive to external perturbation or rotation, which 
can be utilized for quantum parameter estimations [38j |. 
These sub-Planck structures are also found most sensi- 
tive against decoherence in diatomic molecular system of 
hydrogen iodide molecule [3q . In our model, it is worth 
to study the rotational sensitivity due to a infinitesimal 
rovibrational coupling in this physical system. As this 
is not the present goal of this work, we are studying it 
in another context [39]. Fig. [3jg) shows the rotational 
effect in the cat state. The CS is revived at t = 0.5 T rev , 
but shifted in its orbit due to rotation. 

We would like to explain another specific case of the 
phase space structures to reveal the intricacies in a bet- 
ter way. Fig. 2] shows the cat state for higher values of 
j = 81. This cat state is separated in momentum space, 
but exactly superposed in real position space. A sin- 
gle interference ripple in Fig. HJa) has dimension 0.07 
a.u. or 3.7 picomctcr. Although the experimental obser- 
vation of small quantum interference structures is very 
challenging, similar interference ripples in picometer scale 
is recently visualized experimentally for I2 molecule [26j j . 
They observed it in the classical time evolution of a vi- 
brational wave packet at time t = 0.25 T rev + 0.25 T c i 
or t = 0.25 T rev + 0.75 T c /. The same picometer struc- 
tures can be observed for a rovibrating CS of I2 molecule 
at a fixed time (t = 0.25 T rev ) with angular momen- 
tum around j = 81. This structure remains stationary 
in their full crossing pathway in the range 74 < j < 88, 
which keeps the flexibility of experiments. 

Now a natural question can be asked whether there is a 
way to know the orientation (f> corresponding to a j value. 
Here, we provide a numerical estimation of this rotation 
angle. We operate an operator U = e lJ °^, where Jq is the 
angular momentum operator, on the initial wave packet 
and obtain 

n' 

71=0 

= x(y,t). (15) 

Then we compute the overlap | < t)\$(y, t) > | 2 , 
numerically and take the value of <p for the largest over- 
lap ~ 1. We repeat it for various angular momenta and 
find the corresponding angles of rotation. In particular, 
we have chosen the rotation of cat state, i.e., at time 
0.25 T rev . The result is depicted in Fig.Oupto a full pe- 
riod of rotation (2tt radians). Two particular rotations 
(j — 60 and j = 81) are indicated in the figure by dotted 
lines. From the nature of the plot, one can infer that 
the wave packet is more phase sensitive at larger j val- 
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ues. We have checked that these rotation angles are in 
good agreement with the phase space pictures. To be 
more precise, let us take a particular case for j = 81. 
We obtained an angle of rotation <f> — 0.5 737T radians for 
j = 81, but the cat state is exactly vertically situated 
(0.57T radians) with respect to the p = line (Fig. BJb)). 
This is because, the initial cat state (j=0, Fig. [DJc)) is 
not exactly horizontal due to the T rev /T c i ratio, which 
is not exactly an integer and imprinting a small phase 
or rotation equivalent to <f> = — 0.0737T radians. Thus, 
for a perfect vertical cat state, one needs the rotation 
of 0.5737T radians, which is exactly the obtained value in 
our analysis. 

IV. CONCLUSION 

We have showed the time evolution of a CS wave packet 
of rotating Morse potential. We are able to manifest the 
rovibrational coupling of a diatomic molecule through 



this CS dynamics in the experimentally allowed param- 
eter regime. We considered a well behaved and appro- 
priate SU(2) coherent state, including few bound states 
of I2 molecule. This can be generated with an avail- 
able pump beam excited at the 10th energy level and for 
angular momenta j = 60 or 81. The rovibrational cou- 
pling is explored in position space, where the physical 
rotation of the CS is not so transparent. Phase space 
structure clearly showed the rotational influences on the 
cat state as well as the compass-like state. Our study 
deprecates the complication of six-dimensional physics in 
phase space for ro-vibrational dynamics. The 1-D rotat- 
ing Morse potential can well capture the rotational effect 
throughout the time evolution in phase space. Moreover, 
we provide a way to identify the amount of rotation of 
the rovibrating CS at arbitrary time and angular mo- 
mentum, which can be very useful to know the behavior 
of the system at certain extend, without going into the 
phase space description. 
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